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In this article, we introduce a new dielectric function that is different even from the 



generalised RPA, first introduced in x>nd-mat/9810043, in that it pays attention to fluctua- 
tions in the momentum distribution. We argued that this is important when the momentum 
distribution is very nonideal. This new dielectric function reduces to the generalised RPA 
OS . when the number-fluctuations are small, and the generalised RPA in turn reduces to the 

OS ■ traditional RPA dielectric function when the momentum distribution is close to ideal. 

> _ 

Let us write the generalised RPA hamiltonian and try and compute the dielectric function. 



Ho = Y, £k "o( k ) - E w E n °( k + q/ 2 M k - q/ 2 ) W 



k q^O k 

H ext (t) = J2(U ext (qt) + U^-qt)) $> q (k) (2) 

oo 
m . 

here rc. q (k) = c k +q/ / 2 c k-q/2 an d ^o(k) = c k c k Let us now write down the quation of motion for n q (k). 
i^n q (k) = ^e k >[n q (k),n (k')] - ^ ^[n q (k),n (k')]n (k' -q) 

OS . k' qVO k' 

g : + E (EMq «) + t^t(-q*))[nq(k),n q /(k')] (3) 

^ : 

Ch , Now, 

o 
o 



[n q (k),ry(k')] = [c k+q /2 Ck -q/ 2 ' c I'+q72 c k'-q72] 



^ . - 4+q/2 C k'-q72 (5 k-q/2,k'+q72 " C I'+q72 Ck -q/ 2,5 k+q/2,k'-q72 ( 4 ) 

One may approximate this as, 

[nq(k),n q /(k')] = [c k+q/2 e k _ q/2 , cj c , +q , /2 c k /_ q / /2 ] 

= [no(k + q/2) - n (k - q/2)]J k jk ^ qi _ q , (5) 

The next approximation would be to replace the number operator by its c-number expectation value. But 
we shall desist from that for the moment. 

[n q (k),n (k )] = n q (k)(<S k > )k _ q/2 - <V k+q/2 ) (6) 

i gl n i( k ) = (ek-q/2 - £k+q/2)"- q (k) - "^r( n o( k - q/2 - q' ) - n (k + q/2 - q'))n q (k) 

qVo 

+ (C^rt(-qt) + C^rt(qt))(no(k + q/2) - n (k - q/2)) (7) 
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(n q (k)) - C/ ext (-qi)C q (k) + [/ e * xt (qt)£> q (k) (8) 
Let us now ignore fluctuations in the momentum distribution. That is, we are allowed to replace 

Mk')n q (k)) = <n (k')>(n q (k)> (9) 

uj U^ti-^C^-co [/ e * xt (qi)D q (k) = (e k _ q / 2 -e k+q/2 )C/ e;ct (-qt)C q (k) + (e k _ q/2 -e k+q/2 );7 ( r xt (qt) J D q (k) 
- E ^((«o(k - q/2 - q )) - <n (k + q/2 - q')))[/ ext (-qt)<7 q (k) 

- E Y"«"o( k - ^ -«»')>- ^°( k + ^ 2 - q')))^xt(qi)^q(k) 

q'/o 

+ (^ext(-qt) + C/;, 4 (qi))(K(k + q/2)) - (n (k - q/2))) (10) 

w C q (k) = (e k _ q/2 - e k+q/2 )C q (k) - J2 Y"«no(k - q/2 - q )) - (n (k + q/2 - q')»C q (k) 

qVo 

+ ((no(k + q/2))-<no(k-q/2)» (11) 



Since, 



I7 e //(qt) = EWq*) + y (p-q)^ext(qi) (12) 



v c_ (k) - y (no(k ~ q/2)) ~ (no(k + q/2)) (w) 

k q k W - ? k+q/2 + ? k-q/ 2 

If we use the definition of the dielectric function we get precisely the WRONG ANSWER !!! 

, v U ext (qt) 1 
e^iJOiv G (q, co) = — — - = - (no(k _ q/2)) _ (no(k+q/2)) (15) 

/M4; lf V^k w-e k+q / 2 +?k- q /2 

The reason is because we have been very cavalier in our treatment of correlations. Parenthetically we note 
that when I said this is the wrong answer, I meant it is not the RPA dielectric function. But then so what ?? 
Nobody told me RPA was a controlled approximation. But then we must press on, and try to find a clever 
way of reproducing the RPA dielectric function. One has to consider the full Hamiltonian when dealing with 
the dielectric function rather than just H . Let us now try and do this. 

*^<(k) = (ek-q/2 - e k+ q/ 2 K(k) + y ^Kj(k),P q <K_ q , + y ^P q 'K(k),p_ q ,] 

qVo qVo 

+ y (CWq't) + ^ t (-q't))[«q(k),P q '] (16) 
qVo 
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K( k ),Pq'] = H[ C ic+q/2 Ck -q/ 2 ' C I'+q'/2 C k'-q72] 



k' 

= 5Z C k+q/2 C k'-q72 5 k-q/2,k'+q72 ~ X] C k' +q' / 2 C k-q/2<W q / 2: k' -q'/2 
k' k' 

« * q ,- q ' K(k + q/2) - n (k - q/2)] (17) 

*^<( k ) = ( £ k- q /2 - ek+q/ 2 )<(k) + ^Mk + q/2) - n (k - q/2))p q 

+ (Uext(-qt) + f/ e * xt (qi))(n (k + q/2) - n (k - q/2)) (18) 

Let us make a first pass at the computation of the dielectric function. Here, we make use of mean-field 
theory, that is, replace (no(k )p q ) = (no(k ))(p q ) 

w (n* (k)> - (e k _ q/2 - e k+q/2 )«(k)) + ^«n (k + q/2)) - (n (k - q/2)))(p q ) 

+ (^ext(-qi) + C^ rt (qi))«n (k + q/2)) - (n (k - q/2))) (19) 
, _ (n (k + q/2))-(n (k-q/2)) / 

\ n ql K J/ " ,, , \Pq/ 

V W - £k- q /2 + £k+q/2 

+ (t /.„(- qt ) + iCM) ( "° (k + " /2)> ' l""" 1 ' " /2)> (20) 

W — £k-q/2 + e k+q/2 



This means, 



<P- q > = CWq*)%# (2D 



p. x \- (no(k-q/2))-(n (k + q/2)) 

rr V - £ k +q/2 + £ k-q/2 



(22) 



e(q,w) = l-^P (q,w) (23) 

From this and the fact that 

/ \ C^ea;t(qi) , n , n .s 

e g - R p A (ci, u) = - 7— = e(q, w) (24) 

Ueff{qfi) 

Next we would like to include fluctuations. Let us do this differently this time via the use of the BBGKY 
heirarchy. 

4«(k)> = (e k _ q/2 - 6 k+q/2 )(n q (k)) + ^«no(k + q/2)) - <n„(k - q/2)))(p q ) 
+ ^(F 2yl (k + q/2,q) - F 2A (k - q/2, q)) 
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+ (CW-qt) + U: xt (qt))(n a (k + q/2) - n (k - q/2)) (25) 

Here, 

F 2A (k,q;t) = (no(k')pq) - («o(k')>^) (26) 

F 2 (k';k,q;t) = (n (k')n q (k)) - (n (k'))<n q (k)) (27) 

F 2A (k',q;i) -^F 2 (k';k,q;t) (28) 

k 

*-F 2 (k ; k, q; i) = (e k _ q/2 - e k+q/2 )F 2 (k ; k, q; t) + ^ (iV(k , k + q/2) - iV(k' , k - q/2)) (p q ) 

+ ^((n (k + q/2)) - (no(k - q/2)»F 2A (k',q; i) + (I/^-qi) + ?7 e * xt (qt))(7V(k',k + q/2) - JV(k', k - q/2)) 

(29) 

Let us write, 

F 2 (k';k,q;i) = (7 e;Et (-qt)F 2 , a (k'; k, q) + C/ e * rf (qt)F 2 , b (k' ; k, q) (30) 
(p q ) = f/ ext (-qt)(p q ) + f/; xt (qi)(p q ) (31) 

Also define, 

7V(k,k') = (no(k)no(k')) - (n (k))(n (k')) (32) 
cj F 2 , a (k'; k, q) = (e k _ q/2 - e k+q/2 )F 2 , a (k'; k, q) + ^(W(k',k + q/2) - iV(k',k - q/2))(p q ) 

+ ^((no(k + q/2)) - (no(k - q/2)))f 2 a 4 (k',q) + (JV(k',k + q/2) - JV(k',k - q/2)) (33) 

-u,^ 2 , 6 (k^k,q) = (e k _ q/2 -e k+q/2 ^ 

+ ^((n (k + q/2)) - (no(k - q/2)))^(k',q) + (JV(k',k + q/2) - JV(k',k - q/2)) (34) 

I \jpan' s « q ^ iV(k',k + q/2)-^(k , ,k-q/2) / , 

e(q,w) F 2jl (k,q) = ^ — — — (p ) 

V 7 ^ w - e k _ q/2 + £ k +q/2 

+ y iV(k',k + q/2)-iV(k',k-q/2) (gg) 

^ W - £k-q/2 + £k+q/2 

a- C q (k) = (e k _ q/2 -e k+q/2 )C q (k) + ^((n (k+q^ 

+ ((no(k + q/2))-<no(k-q/2)» (36) 



r (v , _ "g (no(k + q/2))-(no(k-q/2)) / ^ 

V U>- £ k -q/2 + Ck+q/2 

(no(k + q/2))-(no(k-q/2)> 



+ - 



e k-q/2 + e k+q/2 



, ^ Fj A (k + q/2,(j) -F 2 a A (k-q/2,q) 

V U — £k-q/2 + £ k+q/2 

After all this, it may be shown that the overall dielectric function including possible fluctuations in the 
momentum distribution is given by, 

£e// (q,c) = e s _Wq,") - P2(q ;"L (38) 

Here, 

^ JV(k + q/2, k' + q/2) - JV(k - q/2, k' + q/2) - jV(k + q/2, k' - q/2) + 7V(k - q/2, k' - q/2) 

^ (W - £ k -q/2 + £k+q/2)(w - <V-q/2 + e k'+q/ 2 ) 

(39) 

/ \ , , "q V- (n (k + q/2))-(n (k-q/2)) 

6 f-« M = 1 + — .L — ^37 , 1 1 - — : — ( 4 °) 



V u 



£ k+q/2 + £ k-q/2 
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